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Abstract. Using a direct approach the return map near a focus 
of a planar vector field with nilpotent linear part is found as a 
convergent power series which is a perturbation of the identity 
and whose terms can be calculated iteratively. The first nontrivial 
coefficient is the value of an Abelian integral, and the following 
ones are explicitly given as iterated integrals. 

1. Introduction 

The study of planar vector fields has been the subject of intense 
research, particularly in connection to Hilbert's 16th Problem. Signifi- 
cant progress has been made in the geometric theory of these fields, as 
well as in bifurcation theory, normal forms, foliations, and the study of 
Abelian integrals [6J, [20] . 

The Poincare first return maps have been studied in view of their 
relevance for establishing the existence of closed orbits, and also due to 
their large number of applications (see for example [12] and references 
therein), and also in connection to o-minimality [IB] . 

A fundamental result concerns the asymptotic form of return maps 
states that if the singular points of a C°° vector field are algebraically 
isolated, there exists a semitransversal arc such that the return map 
admits an asymptotic expansion is positive powers of x and logs (with 
the first term linear), or has its principal part a finite composition of 
powers and exponentials [13], [T7] . 

In the case when the linear part of the vector field has non-zero eigen- 
values there are important results containing the return map P] , [3] , [1] , [3] , 
[Z], [S], il8j. Results are also available for perturbations of Hamiltonians 
[S], [ID] and for perturbations of integrable systems [TT]. On the other 
hand, there are few results available in the general setting [1], [TJ], |15j . 

The present paper studies an example of a field with nilpotent linear 
part, near a focus. The main goal is to establish techniques that allow 
to deduce the return map as a suitable series which can be calculated 
algorithmically and can be used in numerical calculations. 
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2. Main result 
The paper studies the return map for the system 
(1) X = -Y, Y = X^ -Y^ 

which has a nilpotent hnear part (both eigenvalues are zero). This is 
one of the simplest examples of systems in this class [19], and for which 
there are (to the author's knowledge) no methods available to generate 
the return map. 

The main result is the following: 

Proposition 1. Let e with < e < eo small enough. 

The solution of (Q]) satisfying X(0) = e, F (0) = first returns to the 

positive X-axis at the value X satisfying 

oo 

3n+l 



n=l 

which is a convergent series. 

The coefficients X„ can be calculated iteratively. In particular, 

Xi = -2^/2 = 16 Cl^ X3 = -2^/2 (^20 C2C1 + 8 C3 + 49 ci^) 

where Cn = f„(l) with Vn given by 

(2) v,{i) = t t'p{tfi^ dt, V2{i) = -lr^ tt^p{tf/\,{t)dt 
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(3) ^3(0 = e-^ / t 

Jo 



3 ...1/2 , 3 vi^jt) 



dt 



where 

(4) p{t) =A-6t'^ + At^ -t^ 

3. Proof of Proposition[T] 

The proof of Proposition!!] also provides an algorithm for calculating 
iteratively the coefficients 

3.1. Normalization. It is convenient to normalize the variables X, Y, t 
so that the constant e appears as a small parameter in the equation: 
with 

X = tx, Y = T = 2-^^\t 

the system ([1]) becomes 

dX ^ dy n n 

(5) — = -2y, ^ = Ax'- ay^ 



FIRST RETURN MAP 



3 




a = 2^/2,3 



While the initial condition X(0) = e becomes x(0) = 1, it is useful to 
study solutions of ([5]) with the more general initial condition a;(0) = t] 
with 7] in a neighborhood of 1. 

Remark. The system ([5]) has the form dH + auo = with H = 
+ , uj = dx and a a small parameter. The recent result [9J gives 
a formalism for finding the return map for this type of systems; see 
also [llj. The present construction is concrete, explicit, and suitable 
for numerical calculations. 

3.2. General behavior of solutions of 1^. Let a > 0. 

Note the following Lyapunov function for ([5]): 



(7) Lix,y)=y^ + x\ with — L (a;(r), y(r)) = -2« |/(r)^ < 



Since the set {y = 0} contains no trajectories besides the origin 
(which is the only equilibrium point of ([5])), then the origin is asymp- 
totically stable by the Krasovskii-LaSalle principle. 

Consider the solution of ([5]) with the initial condition x(0) = r], y{0) = 
for some 7] G [1/2,3/2]. 

Since y'{0) > and x'(0) = 0, x"(0) < then y increases and x 
decreases for small r > 0. This monotony must change due to ([7]), and 
this can happen only at some point where y = or where Ax^ = ay^, 
whichever comes first. Since y increases, then the first occurrence is a 
point where 4x^ = ay^. At this point x' < so x continues to decrease, 
while y" = —2Ax'^y < so y has a maximum, and will continue by 
decreasing. Again, the monotony must change due to ([7]) and the path 
cannot cross again the line = a^^^y before the monotony of x 

changes, therefore the the next change of monotony happens for y = 0, 
a point where, therefore x < 0; denote this value of x by —fj. 

Denote by Xo(t),?/o(t) the solution for a = 0: Xq = —2y, ?/q = 4xq 
and Xo(0) = r], yo{0) = 0. Therefore x^ + y^ = rj'^. We have x(r) = 
Xo{t) + 0(q;), yij-) = yo^r) + 0{a) therefore fj = r] + 0{a). 

Similar arguments show that the solution x{T),y{T) of ([5]) continues 
to turn around the origin, crossing again the positive x-axis at a point 
fj = 1] + 0{a) (and, of course, f] < r] hj ([7])). 

Solutions (x(r),|/(r)) of provide smooth parametrizations for so- 
lutions y{x) of 



d 

dx 



(y^) = -4x3 + ay^ 
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Note that following the path (x(r), ?/(r)) one full rotation around the 
origin corresponds to considering a positive solution of , followed by 
a negative one. 

3.3. Positive solutions of (IHD for x > 0. LemmaE] shows that there 
exists a unique solution ?/ > of (IHl) so that y{ri) = 0, that this solu- 
tion is defined for x G [0, 77] and establishes an iterative procedure for 
calculating this solution. 

Substituting y = v^l'^ in equation ([8]) we obtain 

(9) ^ = -4x3 + au^l'' 

ax 

Lemma 2. There exists 60 > independent of rj and a so that so that 
the following holds. 

Let 7] > 0. For any a with |a;| < ao < 5ori~^ /2 equation ^ with the 
condition uirj) = has a unique solution u = u{x;a,ri) for x G [O,//]. 
We have u{x; a,ri) > for x E [0,ri) , a > and u{x] a, rj) is analytic 
in a and rj for \a\ < and rj > 0. 

Remark. We will use the results of Lemma[2] only for rj such that 
1/2 < 7] < 3/2. In this case (by lowering oq) we can take of Lemma[2] 
independent of 77 by taking 

(10) ao < 5o/2 min rj'^ = 4727^50/2 = cq^o 
Proof of LemmalE 

Local analysis shows that solutions of ([9]) satisfying u{ri) = have 
an expansion in integer and half-integer powers of r] — x and we have 

1 fi 

(11) u(x) = r]^-x^ ar]^/\7]-xf^l + o(l)) (x ^ r]-) 

5 

which inspires the following substitutions. 
Denote 

(with the usual branch of the square root for x/r/ < 1) and let 

u (x) = r]^ — x^ — 77^ (r] — x) f (^) 

Note that if u{ri) = then necessarily f (0) = by ( fTTjl . 
Equation (jH]) becomes 

(12) ^-^v{0 + 2v (0 = 5 e [piO - V (0 
where 

(13) 5 = 2ri^a 
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and the polynomial p is given by (jll). (Note that p(^) G [1,4] for 

ee[0,i].) 

Lemma[2] follows if we show the following: 

Lemma 3. These exists 6q > so that for any 6 with \6\ < 6q equation 
has a unique solution v = v{^; 6) for G [0, 1] so that v{0) = 0. 
Moreover, v{^; ■) is analytic for 5 G C with \6\ < 6o and the terms of 
its power series 

(14) v{^;6) = J2S''MO 

n>l 

can be calculated recursively; in particular, the first terms are 

To prove Lemma[3] multiply (fT2|) by ^ and integrate; we obtain that 
is a fixed point {v = J'[v]) for the operator 
(15) 

Jv (0 = 5r'/ t^[p{t)-v{t)f'dt = 6e [ s^[p{^s)-v{^s)f'ds 
Jo Jo 

Let B be the Banach space of functions f{C,] 6) continuous for G [0, 1] 
and analytic on the (complex) disk \6\ < 6q, continuous on \6\ < 6o, 
with the norm 

11/11 = sup sup \fi^;6)\ 

«e[0,l] |<5|<<5o 

Let m be a number with < m < 1. Let 5o > be small enough, so 
that 6o < m/y/b and 6o < lO/S/y/b. 

Let Bm be the ball i3„ = {f eB^; \\f\\ < m}. 

We have J' : Bm — > Bm- Indeed, note that for / G Bm we have 

\p{t)-f{t;6)\>\p{t)\-\f{t;6)\>l-m>0 

therefore, since / is analytic in 6, then so is [p{t) — f {t; 6)f'^'^ , and 
therefore so is J'f. Also, if / G Bm then also J'f G Bm because 

\jf (^; 5)\ < \5\ [\p{^s)\ + Ifi^s; 5)\f' ds < '^"^^ < 5oV5 < 

Moreover, the operator JT" is a contraction on Bm- Indeed, using the 
estimate 

\{p-flf'-{p-f2?^'\ < sup \p-f\'/' < 1/1-//^^ + ^^'^' 

we obtain 



2 |/|<m 2 



\jfl-Jf2\<c\\f^-f2\\ withc = 5o^<l 
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Therefore the operator J has a unique fixed point, which is the 
solution t>(,^; 5). 

To obtain the power series (fT^ substitute an expansion v{^\S) = 
^o(0 +En>i ^"^n(0 in (in]). It follows that ^v'q + 2vo = with vo{0) = 
0, therefore fo(0 = 0- 

Substitution of f|T4|) in (p — t>)'^/^ followed by expansion in power 
series in 6 give 

3/2 

(p-^;)3/2=p3 



V n>l ^ / \ n>l / 



where i?„ = i?„[fi, . . . , f„,p]. In particular, 

3 t;i / 3 3 

-Kl - -- — , -K2 - -77 ^ o 

2 p \ 2 p 8 p^ 

From (fT2|) we obtain the recursive system 

(for > 1 and with Rq = 1), with the only solution with f„(0) = 
given recursively by 

Vn{0=r' P{tf' Rn-l{t) dt 

Jo 

In particular, we have ([2]), ([3]). □ 

The following gathers the conclusions of the present section. 

Corollary 4. There exists > so that for any rj e [1/2, 3/2] and a 
with \a\ < ao equation ^ has a unique solution y{x) on [0,1]] satisfying 
y{r]) = and y > on [0, rj) for a > 0. 
Moreover, this solution has the form 

(16) y = (p{x; a,r]) = i]'^ — — r]^ {i] — x) v iJ 1 ■,2ri^a\ 

with V = f (^, 6) a solution of ^^). The map a t— > f (^; 2rfQ.) is analytic 
for \a\ < ao. 

□ 
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3.4. Solutions of (I8j) in other quadrants and matching. 

3.4.1. Solutions in other quadrants. We found an expression for the 
solution y{x) of ([8]) for x > and y > 0. In a similar way, expressions 
in the other quadrants can be found. However, taking advantage of the 
discrete symmetries of equation ([8]), these solutions can be immediately 
written down as follows. 

Let rj G [1/2,3/2], a with |a| < Oq. 

(i) Let yi = (f){x]a,ri) the solution ([TB|) . defined for x G [0,//], with 
yiiv) = ^iid > for a > 0. 

Then: 

(ii) the function 1/2 = (p{—x; —Oi^ri) is also a solution of ([8]), defined 
for X G [— r], 0]; we have y2{—v) = and 1/2 > for a > 0. 

(iii) The function 7/3 = — 0(— x; a,ri) is a solution of ([8]), defined for 
X G [— ry, 0]. We have ysi—rj) = and t/s < for a > 0. 

(iv) The function 1/4 = — —a,ri) is a solution of ([8]), defined for 
X G [0, rj] and we have ^4(77) = and ?/4 < for a > 0. 

3.4.2. Matching at the positive y-axis. Let ri,fj ^ [1/2,3/4] and let 
yi{x) = (l){x;a,ri) solution of ([H]) as in (i), for x G [0,ri] and y2{x) = 
(/)(— x; —a,//) solution as in (ii), for x G [— r/,0]. 

The following Lemma finds fj so that yi{0) = 1/2(0), therefore so that 
yi is the continuation of y2- 

Lemma 5. Let \a\ < cto with ao satisfying /[U^) . Let rj so that \ri — l\ < 
Ci\a\ with ao small enough so that CiOq < 1/2. 
There exists a unique fj = rj + 0{a) so that 

(17) 0(0; a, 77) = 0(0; -a, 77) 

and fj depends analytically on rj and a for \a\ < ai for ai small enough. 

We have |?7 — 1| < C2\a\ for some C2 > and 

(18) 

fl = n- T]^cia + 2 a^ri^ci^ - a^r]^° C3 + y + 10 C2Ci^ + O(a^) 

where Cn = fn(l) with f„,(0 given by 

Proof 
Let 

(19) F(r7,r7,a) = 0(0; -a,//) - 0(O;a,r7) 

which is a function analytic in (f/,?7,a) by Lemma[2]and relation ( fT3|l . 
We have F{r], 77, 0) = and 

dF 

— if], r], 0) = -4r/^ [1 - ?;(1; 0)] = -Ar]' ^ 
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therefore the imphcit equation F{fj,ri,a) = determines fj = fi{ri,a) 
as an analytic function of a, rj for a smaU. 
We have 



I ~ I ^ I I I I 'II 

\ri — ri\ < \a\ sup | — |=C;^|a| 



therefore 

\v — M — \v ~ v\ + \v ~ M ^ (ci + c'i)|a| = C2\(y\ 

The expansion of fj in power series of a is found as follows. Using 
f|T6|) equation f|T7|) becomes 



f]^ - r]\ (l; -2 rf'a) = r]^ - r]^v (l; 2 rfa) 

where, substituting ([H]) and fj = J2n>o ^"Vn followed by power series 
expansion in a we obtain (ITSl) . □ 



3.4.3. Matching at the negative y- axis. Let r], G [1/2, 3/4] and given 
by LemmaO Consider the solution = —(f){—x; a,fj) as in (iii), for 
X e [—fj, 0], with ysi—fj) = 0. Therefore y^ is the continuation of y2- 
Let y4,{x) = — —a, fj) be a solution of ([HD as in (iv), for x G [0, 
The following Lemma finds fj so that 1/3(0) = 1/4(0), therefore so that 
1/4 is the continuation of y^: 

Lemma 6. Let \a\ < ai and \fj — 1\ < C2|a| with ai small enough so 
thatfj e [1/2,3/2]. 

There exists a unique fi > so that 0(0; —a, fj) = 0(0; a, fj) and fj 
depends analytically on fj and a for \a\ < for 0^2 small enough. We 
have 1^ — 1| < cslal for some C3 > and 
(20) 

fj = fj- fcia + 2 a^fj'^ci^ - a^^^ ^3 + y ci^ + 10 csCi^ + O(a^) 

where Cn = fn(l) with Vn{C,) given by 
Proof. 

We need to find fj = fj{fj, a) so that 

F{fj, fj, a) = 0(0; -a, fj) - 0(0; a,fj) = 



Note that the function F above is the same as (fT9l) . By Lemma[5] 
the present Lemma follows. □ 
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3.5. The first return map. Let rj G [1/2, 3/2] and a2 as in LemmaEl 
Tiien fj given by Lemma[6] is the first return to the positive x-axis of 
the solution with x(0) = ri,y{0) = and it is analytic in a and rj, 
therefore, by f|T3|) . it is analytic in e for fixed 1]. 

Combining (fTSl) and ( !20l) we obtain r/ as a convergent power series 
in rj, with coefficients dependent on rj, whose first terms are 

(21) f] = T] — 2 rfc\OL + 8 -if c\c? 

- r/^° (20 C2C1 + 8 C3 + 49 c^) + O (a^) 

To obtain the point X where the solution of ([T]) with X(0) = e > 
and 1^(0) = first returns to the positive X-axis let a = 2^/^e'^ and 
multiply fl^T]) by e (since we have X = ex) and, finally, let r/ = 1. We 

obtain that X is analytic in e for small e and 



1 = e - 2^/2 cie^ + 16 c^h' - 2^/^ (20 C2C1 + 8 C3 + 49 ci=^) 6^° + O {e'') 



where Cn = fn(l) with t>„ given by ([2]), (j3]). 

Remark. The first coefficient of the return map (12T!) is, up to a sign, 
the Melnikov integral of the system ([5]), see[9j; of course, the present 
results are in agreement with this fact (see §l]for details). 

Acknowledgement. The author is grateful to Chris Miller for sug- 
gesting the problem. 



With the notation H = + x'^ and u = dx the Melnikov integral 
of ([5]) is the quantity M(T) = f^^j^u where T > 0. If T is the 
parametrization for the restriction of H to the half-line x > (which 
means, in the notations used in the present paper, that T = r/"^) then 
the return map of ([5]) has the form T \—>- T — aM{T) + 0{a^) [9]. 

For the present system we have 



Taking the fourth root in the return map of T we obtain fl2T!) . 



4. Appendix 
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